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Abstract 

The liquid-gas phase transition in finite nuclei is studied in a heated 
liquid-drop model where the drop is assumed to be in thermodynamic 
equilibrium with the vapour emanated from it. Changing pressure 
along the liquid-gas coexistence line of the systems, symmetric or asym- 
metric, suggests that the phase transition is a continuous one. This 
is further corroborated from the study of the thermal evolution of the 
entropy at constant pressure. 

PACS Number(s): 25.70.-z, 21.65.+f, 24.10.Pa 
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The study of liquid-gas phase transition in finite nuclear systems is of consid- 
erable contemporary interest U,0,[^|§,[6|,[7],|?|[| . Experimental analyses of the accu- 
mulated data on multifragmentation and caloric curves show compelling evidence 
of such a transition. Phase transitions are normally signalled by peaks in the spe- 
cific heat at constant volume Cy with rise in temperature. Theoretical models of 
different genres, such as the microcanonical [|l(J or the canonical [11] description 
of multifragmentation, the lattice-gas model |?|,|8|,|T2"|| or even the microscopic treat- 
ment in a relativistic or a nonrelativistic [[| Thomas- Fermi framework support 
such a structure in the heat capacity. A clear idea about the subtle details of the 
liquid-gas phase transition in finite nuclei, however, has not emerged yet. Confu- 
sion remains about whether the system has evolved dynamically through the critical 
point []6|)[l4|,|15| ; a coherent picture about the order of the phase transition is also 



missing. Analyses of the EOS group [16] and the ISiS group JT^| in the scaling model 
give strong circumstantial evidence for a continuous (second order) phase transition. 
Calculations performed in a mean field model |5j also lead to a similar conclusion. 
Predictions from the lattice-gas models [0,0] are, however, compatible with a first 
order transition. 

Symmetric infinite nuclear matter is effectively a one-component system; with 
heating, it undergoes a first order phase transition. On the other hand, in case 
of two-component asymmetric nuclear matter, as was shown in a comprehensive 



analysis by Midler and Serot [T9|] in a relativistic mean-field framework, the separate 
conservation of the neutron and proton number densities leads the system to a 
continuous phase transition over a finite temperature interval. This is also supported 
in the nonrelativistic calculations by Kolomietz et al |20| . Unlike symmetric nuclear 
matter, even a finite symmetric nucleus (N = Z) behaves like a two-component 
system as the Coulomb interaction lifts the isospin degeneracy. A finite nucleus is 
then expected to undergo a continuous liquid-gas phase transition, if at all. The 
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conflicting predictions from the previous model analyses do not have realistic inputs 
of nuclear physics as relevant for a quantum system of interacting fermions with a 
short-range nuclear force as also with the long-range Coulomb interaction. In this 
communication, we focus on the nature of the phase transition once specific features 
concerning an atomic nucleus are properly taken into account. 

For our study, we choose a representative system, namely, Rhenium with A— 186 
and Z— 75. Such a nucleus is likely to be formed from the reaction 124 Sn + 124 Sn at 
energies around 50 MeV per nucleon GTJ after some nucleons have left out from the 



reaction zone as preequilibrium particles. We also investigate 150 i?e to explore the 
isospin asymmetry effects in phase transition. The nucleus is viewed as a spherical 
liquid drop with asymmetry X defined as (N — Z )/A where A = N + Z , 
is the mass number of the total system. At a finite temperature, we assume the 
depleted nucleus to be enveloped by its own vapour and the system to be in complete 
thermodynamic equilibrium conserving the total number of neutrons and protons. 
The nucleon distributions in the liquid and gas are assumed to be uniform in each 
phase. This definition allows to explore the liquid-gas coexistence region for a finite 
system in close analogy with bulk nuclear matter. 

In absence of a well-defined way to write the energy density functional of a finite 
nucleus in terms of volume, surface, symmetry and Coulomb terms, we write the 
free energy of the nuclear system at temperature T in the single phase as 

F = A f nm (p, X ) + F C + F surf (1) 

where f nm (p, X ) is the free energy per particle of infinite nuclear matter at density 
p with asymmetry Xo at the same temperature T, Fc the Coulomb free energy 
and F sur f the temperature and asymmetry dependent surface free energy. In the 
liquid-gas coexistence region, the free energy is given by 

F co = F l + F 9 , (2) 



where the liquid free energy is, 

F l = A l f nm {p\X l ) + F l c + F l SUTf) (3) 

and the free energy of the emanated gas is 

F 9 = A 9 f nm (P 9 ,X 9 ) + F^. (4) 

Here A 1 and A 9 are the number of nucleons in the liquid and the gas phase, p , p 9 and 
X 1 , X 9 are the corresponding density and asymmetry. The free energy of infinite 
nuclear matter is evaluated in the finite temperature Thomas-Fermi framework with 



a modified Seyler-Blanchard interaction |p2| . The Coulomb free energies for the 



liquid and the gas, Fq and Fq are calculated corresponding to a uniform charged 
sphere and a spherical-shell, respectively. For simplicity, their mutual interaction is 
neglected. The surface free energy of the liquid part F l sur f is taken as |23| 



F l surf = a{X\T){A l f/\ (5) 

where the surface energy coefficient is 

a(X, T) = [a(X = 0) - a s X 2 ] [1 + 1.5T/TJ [1 - T/T c f 2 . (6) 

Here, a(X = 0) = 18 MeV, a s = 28.66 MeV and T c , the critical temperature of the 
symmetric nuclear matter is 15 MeV. The surface energy coefficient decreases with 
density; we neglect its density dependence for the liquid part. Since the gas density 
is very low, its surface energy is neglected. The total surface energy of the liquid is 

surf surf ' ^surf \ ) 

where the total surface entropy is given by 

S Lrf = ~ ( gj 1 ) ■ ( 8 ) 
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The total entropy So of the system is then calculated as 

So = A l s l nm (p\ X 1 ) + A^Jp 9 , X 9 ) + SLf (9) 

The per particle entropy s nm of homogeneous nuclear matter is obtained from the 
standard mean-field prescription. 

The chemical potentials /J n and /J z for neutron and proton in the liquid are given 

by 

i _ 

dF l 

A = w . do) 

The liquid pressure is obtained from 

Similar equations follow for the gas phase. For thermodynamic equlibrium between 
the liquid and the gas, the two chemical potentials and the pressures in both phases 
must be the same, i.e, /J n = p 9 n) p) z = p 9 and P l = P 9 . Along the coexistence 
line, the mass of the liquid drop changes. For a chosen A 1 , the quantities p l , p 9 , 
X 1 , X 9 and A 9 are determined by exploiting the three thermodynamic equlibrium 
conditions and the constraints of baryon number and the total isospin conservation: 

Pn,z = \Pn,z + (! - X )pi,z, 

P X = \p l X l + (1 - \)p 9 X 9 . (12) 

Here p is the average nucleon density and A is the liquid volume fraction. 

The isotherms for the nucleus 186 Re at T = 7, 8, 9 and 10 MeV are displayed in 
Fig.l. For comparison, the isotherm for nuclear matter with asymmetry same as 
that of im Re is also shown at T = 10 MeV. The difference between the isotherms 
for the infinite and the finite system is not insignificant. Though the Coulomb and 



surface have opposing effects on the pressure, the former wins over the latter at this 
temperature. The liquid-gas coexistence lines for 186 Re for the four temperatures 
mentioned are shown by the dotted lines. It is seen that the pressure changes along 



the coexistence line, as seen earlier in the case of asymmetric nuclear matter [19 



the slope of the coexistence lines also increases with temperature. It is further noted 
that at a given temperature and asymmetry, because of the Coulomb effect, the slope 
of the coexistence line of a finite nucleus is more compared to that of asymmetric 
nuclear matter. The variation of pressure along the coexistence line is a pointer 
to a continuous phase transition. Unlike nuclear matter, the coexistence lines do 
not extend from the pure gas phase to the pure liquid phase. At relatively lower 
temperatures, it is found that as the system expands, the size of the liquid drop de- 
pletes and reaches a minimum mass beyond which no thermodynamic equilibrium is 
possible; at higher temperatures, with compression the liquid drop attains a limiting 
mass which decreases with increasing temperature as is evident from the figure. As 
an example, the minimum liquid-drop mass at T = 7 MeV is A 1 = 24 (marked as A 
in the figure); at T = 10 MeV, the limiting liquid-drop mass is A 1 = 130 (marked 
as B). This implies that at lower temperatures, a gas of a finite number of nucleons 
when compressed start nucleating with a minimum mass for the seed in order to 
remain in thermodynamic equilibrium. Similarly, at higher temperatures, for the 
coexisting finite system, the evaporated gas should contain a minimum number of 
nucleons. 

The isospin fractionation along the coexistence line for the nucleus 186 Re at 
T = 8 MeV is shown in Fig.2. The system has proton fraction Y = 0.403 (defined 
as Y = Z/A). As the system prepared in the gaseous phase is compressed, the two- 
phase region is encountered at the point A with the emergence of a minimum liquid 
mass (A 1 = 20) at the point B with a proton fraction Yg larger than Yq. With further 
compression, the gas phase depletes from A to C while the liquid phase grows from 
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B to D attaining the total mass Aq and proton fraction Yq. During compression, the 
proton fractions in both phases decrease, but the total proton fraction Y remains 
fixed, as dictated by the conservation of the total isospin given by Eq.([E|). It is 
evident from the figure that the gas phase is more neutron-rich compared to the total 
system while the liquid phase is comparatively neutron-deficient as also observed 
experimentally This feature becomes more prominent with increasing liquid 

mass. In order to explore the asymmetry effect on isospin fractionation in phase 
transition, the calculated results for the symmetric system 150 Re are also displayed 
in the figure. Contrary to the relatively neutron-rich nucleus 186 Re , here it is found 
that the gas phase is proton-rich. In symmetric nuclei, since \x z is greater than fi n 
because of Coulomb interaction, the separation of the gas phase behaves more like 
that of proton-rich nuclear matter. The occurrence of liquid-gas phase transition in 
symmetric medium-heavy nuclei should then lead to preponderance of proton-rich 
isotopes in energetic heavy-ion collisions and this can be tested in experiments. 
The heat capacity per particle Cy for 186 Re at constant volume (defined as 



Fig. 3 at a representative volume V = lOVo which can be interpreted as a freeze-out 
volume. Here Vq is the normal volume of the nucleus calculated with the radius 
parameter r = 1.16 fm. A very broad bump in Cy with a maximum at T ~ 10 
MeV is seen. The system then corresponds to a liquid part with A 1 around 80, 
the rest of the nucleons being in the gas phase. This is contrary to the results in 
the microscopic mean-field calculations obtained earlier 0,[13| at around the same 
freeze-out volume, where a much sharper peak was observed when the system just 
vaporises completely. At very high temperature, the heat capacity saturates at 1.5 
corresponding to a pure classical gas. 

The thermal evolution of the entropy per particle for the nucleus 186 Re at a con- 
stant pressure P = 0.06 MeV fm -3 is shown by the solid line in figure 4. In contrast 




where E* is the total excitation energy of the system) is displayed in 
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to a first order phase transition where the entropy at constant pressure exhibits a 
discontinuity at a particular temperature (the phase transition temperature), here 
the entropy change is continuous. The noticeable rise in entropy in the temper- 
ature range 8.4 to 10.3 MeV is the manifestation of a liquid-gas phase transition 
in this temperature domain for the chosen pressure. The dashed line in the figure 
corresponds to the entropy evolution for the symmetric nucleus 150 Re. This is quite 
similar to that of the asymmetric nucleus considered, but markedly different from 
that of symmetric nuclear matter [jT^J which shows a sharp discontinuity at the tran- 
sition temperature. The striking similarity between the two results shows that the 
finite size and Coulomb effects are the dominant factors in determining the nature 
of phase transition in finite nuclei. The continuous change of entropy at a constant 
pressure plead in favour of the characterisation of the liquid-gas phase transition in 
atomic nuclei, symmetric or asymmetric, as a continuous one. 

The liquid-gas phase transition in finite nuclei with explicit conservation of the 
baryon number and the total isospin has been investigated in this communication 
in a heated liquid drop model. The peaked structure in the heat capacity, though 
broad, signals the occurrence of a liquid-gas phase transition. From the evolution 
of entropy at constant pressure, one sees that the transition occurs over a range of 
temperatures; this strongly suggests that the liquid-gas phase transition in a finite 
nuclear system is continuous. The simplified assumptions in the model may affect the 
results somewhat quantitatively, but the general qualitative features are expected to 
remain unaltered. The thermodynamic concepts may not be very meaningful when 
the number of particles in one of the phases is very small, still this model serves 
as a window to understand the basic features of liquid-gas phase transition in finite 
systems. 
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Figure Captions 



Fig. 1 The isotherms for the system 186 Re at different temperatures as labelled on 
each curve are shown as full lines. The liquid-gas coexistence lines are shown 
by the dotted lines. The long-dashed line refers to the isotherm at T = 10 MeV 
for nuclear matter with asymmetry same as that of 186 Re. The corresponding 
coexistence line is shown by the filled circles. 

Fig. 2 Evolution of proton fraction (Y) along the liquid-gas coexistence line at T = 
8 MeV for the systems as shown. 

Fig. 3 The specific heat capacity Cy as a function of temperature for 186 Re in a 
freeze-out volume 10Vo- 

Fig. 4 Entropy per particle as a function of temperature at a constant pressure for 
the systems 186 Re and 150 Re. 
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